Persistent currents through a quantum impurity: Protection through integrability 
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We consider an integrable model of a one-dimensional mesoscopic ring with the conduction elec- 
trons coupled by a spin exchange to a magnetic impurity. A symmetry analysis based on a Bethe 
Ansatz solution of the model reveals that the current is insensitive to the presence of the impurity. 
We argue that this is true for any integrable impurity-electron interaction, independent of choice 
of physical parameters or couplings. We propose a simple physical picture of how the persistent 
current gets protected by integrability. 
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The physics of quantum impurities has become an im- 
portant chapter in the evolving theory of strongly corre- 
lated matter. The reasons are several: First, quantum 
impurity problems arguably represent the simplest set- 
tings in which to analyze various aspects of electron cor- 
relations. A case in point is the Kondo effect where a 
magnetic impurity induces an effective electron-electron 
interaction that increases as the energy scale is low- 
ered. Various extensions of the original problem have 
opened up entire new fields of investigations, from the 
modeling of correlated transport in DNA molecules, to 
novel scenarios for non- Fermi liquid behaviors Sec- 
ondly, quantum impurity problems are often tractable by 
exact analytical methods, most prominently the Bethe 
Ansatz technique which exploits the integrability of the 
paradigmatic Kondo- and Anderson models Q ■ This has 
proven immensely useful, with the exact results serving as 
"benchmarks" for more conventional, numerical or per- 
turbative methods. 

Most importantly, progress in experiments on meso- 
scopic and nanoscale systems now enables controlled 
studies of a single quantum impurity interacting with 
conduction electrons. In groundbreaking experiments 
in the late 90s 0], semiconductor quantum dots, con- 
nected capacitively to a gate and via tunnel junctions 
to electrodes, were shown to exhibit a tunable Kondo 
effect: Below a characteristic temperature Tk, a single 
electron occupying the highest spin-degenerate level of 
the dot forms a spin-singlet with electrons in the leads, 
producing a Kondo resonance at the Fermi level. This 
and subsequent developments [4| have turned quantum 
impurity physics into an essential piece also of modern 
nanoscience. 

An interesting problem in this context is how the 
charge persistent current (PC) in a mesoscopic ring cou- 
pled to a quantum dot is affected by a Kondo reso- 
nance. A PC is the equilibrium response to a magnetic 
Aharonov-Bohm flux piercing the ring H . It requires for 
its existence that an electron maintains its phase coher- 
ence while encircling the ring, and is thus expected to be 



sensitive to scattering off the quantum dot. In a previous 
study it was found that the PC is amenable to a Bethe 
Ansatz analysis when the quantum dot is side-coupled to 
the ring. For certain privileged values of the flux, the 
problem was mapped onto the integrable one-dimensional 
(ID) Kondo model with a linearized dispersion. Contrary 
to expectation it was found that the Kondo impurity that 
represents the dot has no effect on the persistent current. 
While this result conforms with those of some other au- 
thors @, HI, a well-controlled RG analysis [§] together 
with large-scale numerics [l(| strongly suggest that the 
PC in fact vanishes when the ring is larger than the 
Kondo screening cloud (other related work includes Refs. 
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I2j). This raised doubts about the applicability of 
a Bethe Ansatz approach (l3j . Having served for many 
years as a work horse in the study of bulk quantum im- 
purity physics, the ID integrable Kondo model was now 
perceived to suffer from a difficulty when applied to this 
particular problem: Its linear dispersion relation, in ad- 
dition to decoupling spin and charge degrees of freedom, 
enforces a non-standard procedure for extracting the PC 
from the finite-size spectrum. It was suggested that these 
features likely explain the failure to obtain an effect from 
the impurity on the PC [l3[ . 

In an attempt to shed light on this intriguing issue, 
we investigate, in this Letter, the influence of a local 
magnetic moment on the PC in a mesoscopic ring, us- 
ing an integrable model with a a non-linear dispersion 
relation for the electrons. As in Ref. the impurity 
is coupled to the ring in such a way that the ring is un- 
affected when the coupling to the impurity is switched 
off. Unlike the analysis in Ref. [6[, however, we do not 
linearize the electronic spectrum, but keep the parabolic 
dispersion of non-relativistic electrons since we want to 
expressly study its possible effect on the PC. Apart from 
the results implied by Refs. the electronic band 

curvature is not known to play any significant role in 
the physics of the Kondo effect. Nevertheless, a small 
number of studies have addressed the question of nonlin- 
earities in the electronic spectrum in the Kondo problem, 
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with a particular eye on how to preserve integrability of 
the model 14 , T^, lj| • In what follows we shall draw on 
some of the insights gained from these studies. 

The basic building blocks that go into the construction 
of an integrable model are the two-particle scattering ma- 
trices Sij [17|. These have to satisfy the Yang-Baxter 
equation 



Sij SikSjk 



SjkSikSij , 



(1) 



the hallmark of integrability [17|. Constructing the 
electron-impurity scattering matrix by the same proce- 
dure as for the ordinary Kondo model 0] but now with 
a quadratic spectrum necessitates for consistency the in- 
troduction of a local potential term in the Hamiltonian: 



V c (x) oc [5'(x+)+5'(x-)]x/\x\ 



(2) 



with x — the location of the impurity. Moreover, 
to satisfy the Yang-Baxter equations (fTJ) for electron- 
impurity and electron-electron scattering, the electrons 
must interact via a local interaction whose strength is 
adapted to the Kondo coupling of the magnetic mo- 
ment. The inclusion of interacting electrons implies a di- 
chotomy: attractive electron-electron interaction neces- 
sitates an antiferromagnetic Kondo coupling, while repul- 
sive electron-electron interaction implies a ferromagnetic 
Kondo coupling. We shall concentrate here on the latter 
case. Since our interest is to study the consequences of 
a non-linear band structure in the framework of an inte- 
grable model, both, the auxiliary potential in eq. ([2]) and 
the dichotomy between electron-electron interaction and 
Kondo coupling, can be easily tolerated. We note in pass- 
ing that a mechanism leading to a ferromagnetic Kondo 
coupling in quantum dots has recently been suggested by 
Silvestrov and Imry [l8j |. 

The first-quantized Hamiltonian on a ring of circum- 
ference L consistent with integrability as outlined above 
is given by [l6j |: 

H = [-dl + (J°i ■ *Q + J')5{xi)} 

i 

i i<j 

where 2 J — — c < and J 1 = — J are required by integra- 
bility. The integrability of the model allows for an exact 
solution, encoded by the Bethe Ansatz equations (BAE) 



Ij/L = z c (kj) 



J 1 j L — Zg(Xj) 
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(4a) 
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Here kj (A 7 ) are rapidities of the charge (spin) degrees 
of freedom, and 9 n (x) = — 2 tan -1 (x/nc). Note that, ex- 
cept for the last term on the right-hand side of eq. (|4bj) 
these are the BAE for the repulsive 5-function Fermi gas 
(lit . This last term, however, encapsulates the contribu- 
tion of the localized magnetic moment. As is well known 
Q , an Aharonov-Bohm flux threading the ring is equiva- 
lent to imposing twisted boundary conditions and in our 
case adds a term proportional to the flux to eq. (|4a|) . We 
shall include such a flux term into our analysis conve- 
niently at a later stage. The quantum numbers Ij and 
J 7 are integers or half-integers depending on the number 
of electrons, N c = N-^+Ni, and the number of down-spin 
electrons, N s — iVj_. Their maximal and minimal values 
are 1^ and J^, respectively. Thus 



N D = J+-J- + 1, 

D c = (/++n/2, 



N s = J+ - J- + 1, (5a) 
D s = (J+ + J")/2. (5b) 



N c and N s constitute the numbers of particles in the 
charge and spin Fermi seas. D c and D s are the numbers 
of electrons and down-spin electrons moved from the left 
to the right Fermi points of their respective Fermi seas. 

It is important to note that in contrast to the standard 
ID Kondo model with a linearized spectrum, charge- and 
spin degrees of freedom are coupled through the BAE 
in ((4}. This suggests that the presence of the magnetic 
impurity may now feed back on the charge sector and 
affect the PC. To find out whether this happens requires 
a careful analysis, to be expounded in what follows. 

Introducing the root density functions p = (p c ,p s ) , 



Pc(k) 



dz c (k) 
dk ' 



dz s {\) 



(6) 



in the charge (c) and spin (s) sectors, we can employ 
the well-known framework [201 ] for extracting the lowest 
order finite-size corrections to the ground-state energy 
in the thermodynamic limit. These are the finite-size 
corrections which determine the PC. Using the Euler- 
Maclaurin formula for converting sums into integrals, one 
can retain finite-size corrections (in principle to arbitrary 
order) when converting the BAE ((4]) into a set of inho- 
mogeneous coupled linear integral equations for the root 
densities [2^]. The phase shifts of eqs. (0} translate, via 
and the Euler-Maclaurin formula, into the integral 
kernel and the inhomogeneity of these integral equations, 
respectively. Given that the phase shifts are odd func- 
tions of their respective arguments, the inhomogencities 
and therefore the solutions of the integral equations also 
attain a certain definite symmetry. Our analysis of the 
PC will eventually rely exclusively on exploiting this sym- 
metry. To achieve a finite-size energy expression correct 
to order 1/L, we introduce the integration limits (fc^A^) 
by 



Lz c (k- 



^±1/2, Lz s (A ± ) = J ± ± 1/2, (7) 
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such that equations ([5]) become 



N c f k+ N s f x+ 

-j- = J dkp c (k), -j- = J dXp s (X), 
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(8b) 
(8c) 
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(9a) 
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dkp c {k)Q 1 / 2 (k) 



x- 



dx Ps {x)e 1 {x) 



(9b) 



The formal solution of the integral equations can then be 
decomposed as 

p(k, A) = poo(k, X) + A) + 

p 1 (fc,A|fc ± ,A ± ) pi(-fc,-A| -fc^-A^) 
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/ ^(/s,A|A: ± ,A ± ) , p 2 (-fc, 



Pc{k- 
-A|-fc* 



-A*) 



p s (A+) 



P-(A") 



(10) 



Here fc ± and A ± play the role of Fermi points of the 
charge and spin excitations. The densities in (fTU|) there- 
fore depend on the numbers and J^, or, via on 
the parameters iV r and D r (r = c, s). The form of (JTUJ) 
depends crucially on the fact that the charge rapidities 
kj enter with odd symmetry into the B AE ((4]) . The term 
Pd/L describes the finite-size contribution of the mag- 
netic moment, with poo, pd, pi and p 2 solving the inte- 
gral equations with appropriate inhomogeneous parts. In 
particular, 
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K 1/2 (X) 



(11) 



with the kernel if 1/2 (A) = ddi/ 2 (X)/dX. 

From these solutions we obtain the ground state energy 
to first order in 1/L, as [20| 
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where Ac = AZ3 C +AL) S is the sum of the finite-size devi- 
ations of D c and D s from their bulk ground state values. 
Eq. IT2|) is obtained by an expansion of the ground state 
energy around the thermodynamic limit L — > 00 which 
in our formalism is tantamount to an expansion in terms 



of (k^ ±A:o) and (X^ ± Ao) around symmetric integration 
limits ko and Ao (cf. (J7J)). This is followed by a transfor- 
mation of variables from ^ and A^ to X r (r = c, s and 
A = AN, AD) evaluated at fc ± = ±fc and A* = ±A , 
thereby incurring as the Jacobian matrix of the trans- 
formation a "dressed charge" matrix [2(| which can be 
shown to obey the same integral equations with the unit 
matrix as inhomogeneity [201 ] . In our case, £ is the func- 
tion that parameterizes this matrix. v c and v s , finally, 
are the Fermi velocities in the charge and spin sectors. 
Note that the leading contribution to the ground-state 
energy due to the magnetic moment is given by edoc- 
This contribution has the character of a boundary term 
2l| . However, the magnetic moment also affects the pa- 
rameters in the 1/L-term, as becomes obvious from the 
decomposition (fTU]) when inserted into © and ©. 

Next, we analyze the finite-size energy to obtain an 
expression for the equilibrium response of the system to 
an externally applied magnetic flux, i.e. the persistent 
current in the presence of the local magnetic moment. 
In fact, the persistent current is precisely determined by 
the finite-size contributions proportional to 1/L in the 
energy (p~2|) . and is obtained by taking the derivative of 
Eq with respect to the external Aharonov-Bohm flux. 
Trading the flux <f> for twisted boundary conditions via a 
gauge transformation leads to an additional shift in the 
number D c of electrons moved from the left to the right 
Fermi points in the charge Fermi sea: AD C — > AD C + <f>. 
Using this replacement in (|12[) . we arrive at the formal 
result for the persistent current: 



ttL 



[AD C + AD S 



(13) 



Now we are in a position to answer the central ques- 
tion of this investigation: How does the presence of the 
magnetic moment, which interacts with the electrons in 
the ring, influence the persistent current? Eq. (|T3|) tells 
us that, to answer this question, we need to analyze the 
effect of the magnetic moment on the parameters AD C 
and AD S . According to © the parts of AD C and AD S 
stemming from the magnetic moment are given by (we 
ignore bulk terms) 



ADt = z d c {0) + ~ 



AD1 



dk- 



dX 



dk 



dX 



pj(k) (14a) 
p d s (X), (14b) 



where the density functions p d {k) and pf(A) are solutions 
of the integral equations with the inhomogeneity pff (cf. 
(jlip ) and we have symmetric integration limits k^ = 
±fco and A^ = ±Ao in the ground state according to our 
discussion after eq. (fT2|) . 

There is, however, no need to explicitly solve the inte- 
gral equations to obtain further insight into the quantities 
AD d and AD d . They follow simply from considering the 
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symmetry of the functions involved. The symmetry prop- 
erties of all functions derived from the integral equations 
follow from the basic odd symmetry of the BA charge ra- 
pidities k and the symmetry of the inhomogcneity. The 
inhomogeneity pdo (cf- (fill) ) i s even, as are all integral 
kernels. Further scrutiny therefore reveals that the sym- 
metries are such that AZ?J? and AZ?f both vanish. E.g. pd 
is an even function in both variables k and A and hence, 
from Q, -Zc(O) = 2g(0) = such that, moreover, from 
(fli]). ADc = and AL>f = 0. Hence there is no influence 
of the magnetic moment on the persistent current. 

We reiterate that our result follows immediately from a 
symmetry analysis of the rapidities and the integral ker- 
nels implied by the BAE ([4]). Since these are generic, the 
result carries over to any model of a quantum impurity 
coupled to electrons with a dispersion relation that is an 
even function of momentum, i.e. for example a parabolic 
(non-relativistic) band, or, if defined on a lattice, a tight- 
binding band. The reason for this universal behaviour of 
integrable quantum impurities is that the details of the 
model do not affect the generic symmetry of the Bethe 
ansatz, as demonstrated by our analysis above. This is 
also consistent with results obtained for the supersym- 
metric t-J model, where the finite-size (oc l/L) contri- 
bution to the energy due to twisted boundary conditions 
was found to be independent of the presence of an inte- 
grable impurity [22 1. 

What is the physics behind this remarkable phe- 
nomenon? We propose that an answer may be con- 
structed as follows: As is well-known, integrable quan- 
tum dynamics in one dimension supports only forward 
scattering [T3|. It is also known that a forward scatter- 
ing phase shift of a free electron wave function incurred 
from a local static potential has no effect on a persistent 
current: As was shown by Gogolin and Prokof'ev [23} 
there is a subtle cancellation (to 0(1/L)) of contribu- 
tions to the persistent current from phase shifted states, 
leading to an expression for the current in terms of the 
Fermi level transition amplitude only. Provided that the 
effect of a quantum impurity on the conduction electrons 
can be faithfully encoded by a potential scatterer (much 
as in Noziere's local Fermi liquid theory of the ordinary 
Kondo effect 24[) - and that this property is not cor- 
rupted on a mesoscopic scale - our result would get an 
elegant and transparent explanation. 

From an experimental point of view one may be con- 
cerned that the protection of the persistent current is not 
robust against any deviation from integrability. That is, 
any small perturbation would make a difference, regard- 
less of the perturbation being relevant or irrelevant in the 
sense of renormalization group theory. Thus, to detect a 
pure protected current will require "fine tuning" of the 
experimental setup so as to make sure that the dynamics 
remains integrable. 

In conclusion, we have demonstrated on quite general 
grounds that there is no influence from a quantum impu- 



rity on the persistent current in a mesoscopic ring when 
the electron-impurity interaction is integrable. We con- 
jecture that this result can be traced back to a cancella- 
tion of phase shifted contributions to the persistent cur- 
rent, in analogy with the simple case of non-interacting 
electrons in the presence of a single forward scattering 
local potential. To put this conjecture on a firm ground, 
and to extract implications for other Aharonov-Bohm (or 
Aharonov-Casher geometries, is an interesting and 
challenging problem. 
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